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1 .  INTRODUCTION 

Work  performed  on  this  contract  has  been  concerned  with 
acoustic  disturbances  which  are  ..elevant  to  the  problem  of  long-range 
detection  and  diagnosis  of  atmospheric  nuclear  explosions.  Our 
earlie:  work  concentrated  on  acoustic  signals  which  originate  in  the 
detonation  bias  *  wave  and  which  can  be  detected  on  the  ground  by 
r  -  '^barographs  or  seismometer  .  More  recently  we  have  investigated 
’}  rising  firet'll  as  a  source  of  acoustic  signals  and  have  turned 
ur  attention  to  acoustic  signals  in  the  ionosphere  which  can  be 
detected  by  electron) -gnetic  methods. 

In  Section  2,  conclusions  reached  during  our  earlier  work 

on  the  contract  are  reviewed.  This  work  is  described  in  detail  in 

* 

our  Semi-Annual  Report.  Three  different  acoustic  phenomena  relevant 
to  nuclear  test  detection  and  diagnostics  were  discussed  in  this 
report:  (1)  The  generation  of  long  period  (3-10  minute)  acoustic- 

gravity  waves  was  analyzed  in  terms  of  a  Lamb  mode  propagation  theory. 
One  result  of  the  analysis  was  that  yield-amplitude  proportionality 
was  predicted  to  break  down  for  very  large  yield  detonations.  (2)  The 
generation  and  propagation  of  short  period  (1-60  seconds)  acoustic 
pulses  was  treated  by  means  of  weak  shock  theory.  Yield  and  height  of 
burst  scaling  laws  were  derived  for  the  far-field  period.  The  depen¬ 
dence  of  period  on  atmospheric  conditions  and  propagation  path  was 
also  discussed.  (3)  The  variation  of  long  range  Rayleigh  wave  ampli¬ 
tude  with  yield  and  height  of  burst  was  discussed  for  atmospheric, 
explosions. 

In  Section  3,  the  relationship  between  hydrodynamic  motions 
caused  by  a  low  altitude  nuclear  explosion  and  subsequent  ionospheric 
disturbances  is  reviewed.  Both  the  upward  going  shock  and  the  rising 
fireball  are  considered  as  hydrodynamic  sources.  It  is  shown  that 
different  portions  of  the  shock  front  may  be  classified  in  terms  of 


B.  L.  Murphy  and  A.  D.  Zalay ,  "Study  of  Acoustic-Gravity  Wave 
Generation  by  Nuclear  Detonations,"  Mt.  Auburn  Research  Associates, 
Inc.  Report  under  AFOSR  Contract  No.  F44620-71-C-0086,  30  September 
1971. 
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the  ionospheric  disturbance  they  create.  The  portion  of  the  shcck 
front  reflected  from  the  100-120  km  altitude  level  produces  distur¬ 
bance  periods  the  order  of  a  minute  for  a  megaton  detonation.  The 
portion  of  the  shock  front  v.hlch  propagates  above  the  100-120  km 
level  is  responsible,  through  a  complex  nonlinear  process,  for 
disturbance  periods  in  excess  of  10  minutes.  It  is  shown  that  the 
fireball  is  most  efficit'\c  in  generating  acoustic-gravity  waves  when 
it  reaches  its  stabilization  altitude  and  approaches  hydrodynamic 
equilibrium  with  the  atmosphere.  Theory  and  laboratory  experiment 
indicate  that  the  fireball  then  produces  a  wave  spectrum  peaked  at 
periods  slightly  longer  than  the  ambient  Brunt-VMs’dia  period  5 
minutes) . 

In  Section  4,  laboratory  experiments  are  described  which 
are  intended  to  simulate  gravity  waves  generated  by  the  rise 
and  stabilization  of  a  fireball.  Observations  were  made  of 
the  internal  waves  generated  in  a  density  scratified  fluid  by  the 
motion  of  a  buoyantly  rising  miscible  fluid.  The  gravity  wave  motions 
were  analyzed  by  observing  the  trajectories  followed  by  small  neutrally 
buoyant  beads  dispersed  throughout  the  stratified  medium  and  good 
agreement  was  found  with  theoretical  predictions.  It  was  confirmed 
that  a  source  of  frequency  u  oscillating  in  a  medium  whose  Brunt- 
Vais'dia  frequency  is  u)q  tends  to  radiate  mainly  at  an  inclination  0 
to  the  vertical  given  by  0  ~  cos  was  °^serve(^  that  the 

rise  time  of  a  buoyant  fluid  to  its  stabilization  height  was  about 
0.85  times  the  Brunt  period,  and  that  the  fluid  typically  undergoes  about 
two  oscillations  before  coming  to  rest  vertically.  The  gravity  wave 
field  for  a  buoyant  release  experiment  was  observed  to  be  inhomogeneous, 
covering  a  range  of  frequencies,  but  in  general  t*>o  >  u>  >  u>  /2.  It  was 
clear  that  gravity  waves  were  generated  mainly  by  the  oscillatory 
motion  of  the  fluid  around  its  stabilization  height,  rather  than  during 
its  rise.  Finally,  it.  was  shown  that  if  measurements  are  made  of  only  two 
characteristics  of  fluid  particle  trajectories,  viz.:  the  inclination 
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of  the  major  axis  of  particle  motion,  and  the  ratio  of  major  to 
minor  axis  for  its  (in  general)  elliptical  orbit,  then  it  is  possible 
to  identify  the  location  of  an  unseen  source  of  internal  gravity  waves. 
Overall  conclusions  reached  are  presented  in  Section  5. 
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2.  REVIEW  OF  PREVIOUS  WORK  ON  THE  CONTRACT 


Three  different  kinds  of  acoustic  phenomena  related  to 
nuclear  test  detection  and  diagnostics  were  discussed  in  the  Serni- 
Annual  Report.  They  are:  (1)  Long  period  (3-iO  minutes)  acoustic- 
gravity  waves  which  propagate  in  the  real  atmosphere's  counterpart  to 
Lamb's  edge  mode.  (2)  High  frequency  (1-60  second  period)  acoustic 
pulses  which  result  from  the  degeneration  of  the  detonation  blast 
wave  into  a  weak  shock.  (3)  Rayleigh  or  surface  seismic  waves  which 
result  when  the  detonation  shock  wave  strikes  the  earth.  Conclusions 
reached  for  each  of  these  phenomena  are  given  below. 

2.1  LAMB  MODE  ANALYSIS 

(a)  The  shock  front  extending  from  the  ground  to  many 
scale  heights  above  the  ground  was  found  to  be  of  almost  equal 
importance  in  determining  the  far-field  Lamb  mode  amplitude.  Because 
of  this,  the  far-field  amplitude  should  be  relatively  independent  of 
phenomena  peculiar  to  any  small  portion  of  the  shock  front.  In 
particular,  ground  effects  such  as  the  presence  of  a  precursor  should 
not  play  a  very  important  role. 

(b)  For  large  yield,  low  altitude  detonations  we  found  a 

f  /2 

far-field  amplitude  dependence  which  varies  as  Yll  +  — j 

where  typically  Y*  is  the  order  of  10  MT  for  a  sea  level  detonation. 

This  dependence  arises  from  a  detailed  consideration  of  the  mechanism 
by  which  the  low  frequency  components  of  the  near-field  (shock)  distur¬ 
bance  become  spatially  separated  from  the  high  frequency  components. 

The  precise  value  of  the  critical  yield  Y*  where  yield-amplitude  pro¬ 
portionality  breaks  down  is  determined  by  atmospheric  conditions  near 
the  burst  as  well  as  the  direction  of  propagation  and  must  be  determined 
by  detailed  numerical  calculations. 

* 

B.  L.  Murphy  and  A.  D.  Zalay,  Study  of  Acoustic-Gravity  Wave 
Generation  by  Nuclear  Detonations,"  Mt.  Auburn  Research  Associates, 
Inc.  Report  under  AFOSR  Contract  No.  F44620-7I-C-0086,  30  September 
1971. 
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2.2  HIGH  FREQUENCY  ACOUSTIC  PERIODS 

(a)  At  large  distances  the  dominant  acoustic  pulse  period 
for  a  1  KT  detonation  was  found  to  be  the  order  of  10  seconds. 

(b)  For  small  yields  the  period  was  found  to  scale  as 
1/3 

(yield)  .  For  large  yields  the  scaling  is  uncertain.  The  inter¬ 
action  of  the  rising  fireball  from  high-yield  detonations  with  the 
upward  going  shock  needs  to  be  investigated  before  the  large  yield 
scaling  can  be  determined. 

(c)  The  range  dependence  of  far-field  periods  was  found 
to  be  extremely  weak,  as  an  absolute  maximum  about  a  17%  increase  in 
going  from  2,000  to  10,000  km.  Pulse  splitting  by  wind  ducts  and  the 
effects  of  other  small  scale  meteorological  phenomena  are  expected  to 
make  the  actual  dependence  even  weaker  than  the  above  value. 

(d)  The  period  dependence  on  height  of  burst  was  also  found 
to  be  extremely  weak.  The  dependence  is  given  approximately  as: 

Ta^l  +  »  where  T  is  the  period  and  z  is  the  burst  height  in  km. 

(e)  Probably  the  most  significant  variation  in  far-field 
period,  other  than  that  caused  by  variations  in  yield,  was  found  to 

be  caused  by  the  direction  of  signal  propagation  relative  to  the  winds 

at  about  50  km  altitude.  For  downwind  propagation  the  signal  reflection 

occurs  at  a  lower  altitude  than  for  crosswind  or  upwind  propagation. 

The  dependence  of  period  on  reflection  altitude  was  approximately  given 
zj28 

by:  T  a  e  ,  where  z  is  the  reflection  altitude  in  kilometers. 

J  o 

A  10  km  variation  in  zq  produces  about  a  40%  variation  in  T.  The 
altitude  zq  can  be  determined  by  ray  tracing  for  the  appropriate 
atmospheric  conditions  and  direction  of  propagation. 

2.3  RAYLEIGH  WAVE  AMPLITUDES 

Our  detailed  conclusions  are  contained  in  the  Semi-Annual 
Report  in  the  form  of  plots  of  source  strength  for  Rayleigh  wave 
excitation  vs  yield  and  height  of  burst.  Some  general  conclusions 
were  as  follows: 
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(a)  The  source  strength,  which  contains  all  the  yield  and 

height  of  burst  dependence  for  the  far-field  amplitudes,  was  found  to 

increase  less  rapidly  with  yield  Y  the  higher  the  burst  altitude.  For 

example,  for  a  Rayleigh  wave  period  of  20  sec  and  for  yields  between 

93 

100  and  500  KT  we  found  S  o  Y*  for  detonations  at  2km  altitude, 

S  a  at  5  km,  S  a  Y*  ^  at  10  km,  and  S  a  at  40  km. 

(b)  For  large  yields  and  short  Rayleigh  wave  periods  t 
(for  example,  Y  >  500  KT  for  t  ■  10  sec)  the  source  strength  was 
found  to  be  a  monotonic  decreasing  function  of  height  of  burst.  For 
larger  yields  or  longer  Rayleigh  wave  periods  S  was  maximized  for  a 
particular  height  of  burst.  The  altitude  at  which  the  maximum  occurred 
increased  with  increasing  x  and  decreasing  Y.  For  Y  =  1  MT  and  x  =  20 
sec  the  maximum  occurred  at  about  45  km  altitude. 

(c)  The  variation  of  S  with  height  of  burst  was  appreciable. 
For  example,  for  x  =  20  sec  and  Y  *  100  KT,  S  increased  by  a  factor 

of  2.3  as  the  height  of  burst  was  raised  from  4  km  to  20  km. 


* 
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3.  MODELING  OF  NUCLEAR  SOURCES  OF  ACOUSTIC-GRAVITY  WAVES* 

3.1  INTRODUCTION 

It  is  known  that  large  yield  nuclear  explosions  result  in 
ionospheric  disturbances  which  can  be  observed  by  electromagnetic  means 
at  great  distances  from  the  explosion  site.  It  is  also  widely  recognized 
that  these  disturbances  are  a  result  of  acoustic-gravity  waves  whose  pro¬ 
pagation  is  supported  by  the  neutral  fluid  component  of  the  ionosphere. 

A  recent  bibliography  contains  many  references  to  observations  and  theory 
relating  to  detonation  produced  disturbances  (THOMAS,  J.  E.,  et  al.,  1972). 

Much  of  the  previous  theoretical  development  deals  with  wave 
propagation  rather  than  source  modeling.  It  is  to  the  latter  topic  that 
this  paper  is  devoted.  In  the  subsequent  discussion  we  outline  how  the 
hydrodynamic  motions  near  a  nuclear  detonation  result  in  the  generation 
of  acoustic-gravity  waves  and  how  the  properties  of  these  waves  are  deter¬ 
mined  by  the  explosion  parameters.  That  is,  how  onset  time,  amplitude, 
and  period  of  the  ionospheric  disturbance  are  related  to  explosion  yield 
and  height  of  burst. 

Because  we  concentrate  on  source  modeling,  we  will  not  discuss 
wave  propagation  except  insofar  as  it  is  necessary  to  do  so  to  relate 
theoretical  predictions  to  observations.  Similarly,  we  do  not  consider 
how  tSie  neutral  fluid  motions  determine  ionospheric  electron  densities 
and  hence  the  ultimate  electromagnetic  effect.  This  difficult  problem 
has  been  treated  by  a  number  of  authors  (for  example,  HOOKE,  W.  H.,  1968). 

Our  treatment  is  limited  to  low  altitude  detonations  where  the 
initial  energy  deposition  may  be  assumed  to  be  local  and  spherically 
symmetric.  For  this  type  of  detonation  both  the  shock  and  the  fireball 
can  generate  acoustic-gravity  waves.  The  mechanisms  which  will  be  dis¬ 
cussed  are  as  follows: 

a.  The  direct  effects  of  the  shock  are  manifest  as  ionospheric 
signals  with  periods  the  order  of  a  minute  for  a  megaton 
detonation,  the  period  being  determined  by  the  positive 
phase  duration  of  the  shock  when  it  reaches  the  ionospheric 

Paper  to  be  presented  at  the  EPP-AGARD  Specialist  Meeting  on  "Effects 
of  Atmospheric  Gravity  Waves  on  Electromagnetic  Wave  Propagation," 
at  Wiesbaden,  Germany  17-21  April  1972. 
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level.  The  acoustic  signal  in  this  case  undergoes  multiple 
reflections  between  the  ground  and  the  base  of  the  ionosphere. 
The  ducting  is  imperfect  at  the  upper  level  and  energy 
continuously  leaks  into  the  ionosphere. 

b.  In  addition,  the  short  period  (order  of  a  minute)  hydro¬ 
dynamic  motions  that  characterize  the  shock  as  it  sntsrs 
the  ionosphere  above  the  burst  evolve  into  long  period 
(10  minutes  and  more)  disturbances  through  a  highly  non¬ 
linear  process  occurring  at  the  100-120  km  altitude 
level.  This  phenomenon,  which  has  been  demonstrated  by 
numerical  calculations  (GREENE,  J.  S.,  Jr.  and  W.  A. 

WHITAKER,  1968)  ,  involves  a  refraction  and  folding  over 
of  the  shock  due  to  the  rapid  change  in  ambient  tem¬ 
perature  at  this  level.  Long  period  wave  propagation  due 
to  this  mechanism  commences  at  10  minutes  or  more  after 
the  detonation.  This  time,  which  is  when  the  effective 
excitation  appears  to  have  begun  according  to  a  distant 
observer,  i3  the  time  required  for  the  shock  to  reach 

the  ionospheric  level  and  for  the  long  period  fluid  motions 
to  become  established. 

c.  The  fireball  is  most  efficient  in  generating  acoustic- 
gravity  waves  when  it  has  reached  its  terminal  or 
stabilization  altitude.  It  then  approaches  hydrodynamic 
equilibrium  with  the  atmosphere  on  a  time  scale  com¬ 
parable  with  the  Brunt-VMisald  period,  the  natural 
oscillation  period  of  the  acmosphere.  This  is  also 
approximately  the  time  required  for  the  fireball  tc 
reach  its  terminal  altitude  and  wave  generation  to 
begin.  For,  as  will  be  shown,  the  fireball  is  not  an 
efficient  emitter  of  waves  during  the  early  portions 

of  its  rise. 

Acoustic-gravity  wave  generation  during  fireball  terminal 
phase  has  been  modeled  in  the  laboratory  and  these  exper¬ 
imental  studies  will  be  described. 


As  shown  below,  for  megaton  range  detonations  the  energy 
available  for  wave  generation  by  the  fireball  is  compar¬ 
able  to  the  energy  available  for  long  period  wave  genera¬ 
tion  by  the  shock.  However,  the  dominant  frequency  of 
the  waves  generated  by  fireball  stabilization  is  only 
slightly  less  than  the  Brunt-VaisSlh  frequency  and  waves 
of  this  frequency  propagate  at  relatively  steep  angles 
to  the  horizontal.  For  this  reason  the  fireball  mechanism 
may  only  be  capable  of  producing  ionospheric  disturbances 
in  the  immediate  vicinity  of  the  burst  location. 

These,  general  remarks  are  discussed  in  more  detail  in  the  succeeding 
sections.  We  first  discuss  the  shock  and  then  the  fireball  mechanism. 

We  then  assess  the  importance  of  these  mechanisms  in  terms  of  the  result¬ 
ing  ionospheric  effects. 

3.2  SHOCK  WAVE  GENERATION  OF  ACOUSTIC-GRAVITY  WAVES 

It  is  useful  to  first  consider  the  propagation  of  the  shock  front 
as  it  is  refracted  by  the  atmosphere.  Different  portions  of  the  shock  are 
refracted  differently.  This  leads  to  a  classification  of  different  portions 
of  the  shock  front  in  terms  of  the  types  of  ionospheric  disturbance  that 
they  create.  ^ 

At  some  distance  from  the  burst  point  the  shock  becomes  weak, 
that  is,  the  relative  overpressure  AP/P  becomes  much  less  than  unity. 

AP  is  the  overpressure  and  P  is  the  ambient  pressure.  For 
example,  for  a  megaton  detonation  AP/P  *  ,1  at  about  9  km  (LEHTO,  D.  L. 
and  R.  A.  LARSON,  1969).  When  the  shock  becomes  weak  the  effects  or 
atmospheric  stratification  and  winds  become  important. 

An  example  of  ray  tracing  calculations  for  a  realistic  (but 
windless)  atmosphere  is  shown  in  Figure  1  (BARRY,  G.,  1963)  for  a  source 
located  at  the  ground.  The  figure  does  not  include  any  nonacoustic  pro¬ 
pagation  effects  such  as  those  arising  near  the  source  where  the  shock 
is  strong.  The  initial  ray  angles  in  Figure  1  start  at  zero  degrees  to 
the  horizontal  and  increase  in  5°  increments.  The  ultimate  behavior  of 
different  portions  of  the  shock  front  is  determined  by  the  initial  angle 
of  propagation  relative  to  the  horizontal  as  indicated  in  Figure  1.  As 
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this  angle  is  increased  the  different  types  of  disturbance  are: 

a.  Acoustic  disturbances  reflected  back  toward  the  earth 
before  reaching  the  ionosphere.  In  Figure  1  this 
corresponds  to  initial  angles  of  less  than  20°.  This 
type  of  signal  may  undergo  repeated  rejections  between 
the  ground  and  some  higher  level,  in  Figure  1  the  40  km 
region  of  the  stratosphere,  where  suitable  temperature 
and  wind  conditions  exist.  When  this  portion  of  the 
shock  front  becomes  weak  it  may  be  repeatedly  split  by 
winds  resulting  in  a  mutipulse  type  of  infrasonic  signal 
(MEECHAM,  W.  C.,  1968). 

b.  The  direct  shock  which  reaches  the  100-120  km  level  of  the 
ionosphere  before  being  reflected  at  the  region  of  very 
sharp  temperature  increase.  This  signal,  which  corresponds 
to  the  interval  from  20°  to  perhaps  35°  in  Figure  1,  can 

be  ducted  between  the  base  of  the  ionosphere  and  the  ground 
as  shown.  As  discussed  below,  leakage  from  the  upper  portion 
of  the  duct  to  the  ionosphere  is  probably  responsible  for 
acous tic-wave  disturbances  observed  in  the  ionosphere 
(BAKER,  D.  M.  and  K .  Davies.  1968  ). 

« 

c.  A  portion  of  the  shock  front  which  propagates  into  the 
thermosphere  beyond  the  100-120  km  level.  This  corresponds 
to  angles  of  more  than  about  35°  In  Figure  1.  Numerical 
calculations  of  the  nonlinear  hydrodynamics  done  by  Greene 
and  Whitaker  (1968),  which  are  discussed  below,  show  that 
this  portion  of  the  shock  loses  much  of  its  energy  in  the 
creation  of  a  horizontally  propagating  disturbance  behind 
the  front. 

This  disturbance  is  believed  to  be  responsible  for  the  very 
long  period  ionospheric  perturbations  observed  at  great 
distances  following  large  yield  explosions. 

Our  concern  in  the  following  is  with  the  portion  of  the  shock 
front  which  reaches  the  100-120  km  altitude  level.  Throughout,  only 
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the  primary  shock  wave  is  considered  and  secondary  shocks  due  to 
reflection  from  the  fireball  or  the  ground  are  neglected. 

We  first  determine  the  properties  of  the  upward  going 
shock  at  the  100-120  km  altitude  level.  We  then  discuss  the  relation¬ 
ship  of  these  shock  properties  to  the  properties  of  short  period  iono¬ 
spheric  disturbances.  Finally  we  discuss  the  nonlinear  process  at  the 
100-120  km  level  which  results  in  the  generation  of  long  period  iono¬ 
spheric  disturbances. 

3.2.1  The  Upward  Going  Shock 

In  this  section  we  will  calculate  the  relative  overpressure 
and  the  positive  phase  duration  of  the  shock  at  the  base  of  the  ionosphere 
as  functions  of  yield  and  height  of  burst.  These  results  are  subsequently 
used  to  relate  ionospheric  disturbance  properties  to  explosion  parameters. 
To  determine  the  properties  of  the  upward  going  shock  a  simplified  scaling 
law  known  as  modified  Sachs  scaling  is  used  (LUTZKY,  M.  and  D.  L.  Lehto, 
1968)  together  with  weak  shock  theory  (REED,  S.  G.,  Jr.,  1959). 

According  to  modified  Sachs  scaling  the  relative  overpressure, 
in  an  inhomogeneous  atmosphere,  at  a  distance  r  from  an  explosion  where 
the  ambient  pressure  is  P(r),  is  just  the  same  as  if  the  explosion  had 
occurred  in  a  homogeneous  atmosphere  with  ambient  pressure  P(r).  That 
is: 


where  Y  is  the  explosion  yield  and  f  is  the  function  which  gives  the 
relative  overpressure  in  a  homogeneous  atmosphere .  An  example  of  the 
application  of  modified  Sachs  scaling  is  shown  in  Figure  2  where  we 
compare  the  predictions  of  this  scaling  with  the  results  of  SAP  and 
SHELL  calculations  for  relative  overpressure  versus  distance.  SAP  is 
a  one  dimensional  Lagrangian  hydrodynamic  code  and  SHELL  is  a  two 
dimensional  Eulerian  hydrodynamic  code.  These  calculations  are  for 
shock  propagation  at  45°  from  the  horizontal  due  to  a  4  MT  isothermal 
sphere  at  5  km  altitude  (GREENE,  J.  S.,  Jr.,  personal  communication, 
1971).  The  modified  Sachs  scaling  curve  has  been  constructed  using 


the  1962  Standard  atmosphere  for  the  ambient  pressure  P(r)  and  the 
homogeneous  atmosphere,  real  air,  calculations  of  Lehto  and  Larson 
(1969)  for  the  function  f  accurring  in  Eq.  (1). 

Notice  that  the  SHELL  calculation  consistently  gives  some¬ 
what  higher  relative  overpressures  below  100  km  altitude  than  does 
SAP.  This  is  because  it  includes  the  effects  of  fireball  rise  which 
SAP  does  not.  The  rising  fireball,  particularly  for  large  yields, 
prevents  the  shock  from  relieving  backwards.  At  very  small  angles 
from  the  horizontal  the  differences  between  SHELL  and  SAP  are  negligible. 
Both  SHELL  and  SAP  show  a  drop  in  relative  overpressure  commencing  at 
about  100  km  altitude.  The  reason  for  this  and  its  relation  to  acoustic- 
gravity  wave  generation  at  the  ionospheric  level  will  be  discussed  short¬ 
ly.  First  we  continue  the  discussion  of  modified  Sachs  scaling  and  the 
upward  shock  propagation. 

It  is  not  entirely  clear  why  modified  Sachs  scaling  works  as 
well  as  it  does.  However,  in  this  example,  as  well  as  in  others  not 
shown  here,  modified  Sachs  scaling  is  definitely  more  accurate  than 
scaling  of  distance  with  respect  to  ambient  pressure  at  the  detonation 
altitude.  Because  modified  Sachs  scaling  is  expected  to  become  more 
inaccurate  with  increasing  distance,  we  use  it  only  to  obtain  starting 
values  near  the  burst  point  for  the  application  of  weak  shock  theory-. 

In  using  iSeak  shock  theory  to  obtain  scaling  laws  for  the  shock  proper¬ 
ties  on  entering  the  ionospheric  level  we  make  two  simplifying  assump¬ 
tions:  First,  that  the  atmosphere  can  be  approximated  as  exponential. 
Second,  that  in  spite  of  refractive  effects,  each  portion  of  the  shock 
front  propagates  independently  with  its  own  effective  radius  of  curva¬ 
ture.  These  assumptions  are  justified  by  the  fact  that  we  treat  shock 
propagation  in  the  atmosphere  below  120  km  altitude  where  the  scale 
height  is  approximately  constant  and  by  the  fact  that  we  only  consider 
the  portion  of  the  shock  front  which  is  not  seriously  refracted  (i.e. 
reflected)  before  reaching  the  altitude  in  question. 

Fo1-  reasons  which  appear  below  we  begin  using  weak  shock  theory 
at  a  distant  ;  where  AP/P  =  APq/Po  first  equals  0.2.  Assuming  that  the 
atmosphere  near  the  burst  may  be  approximated  as  exponential  this  value 
of  the  function  f  occurring  in  Eq.  (1)  corresponds  to  (LEHTO,  D.  L.  and 


R.  A.  LARSON,  1969): 
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where  Y  is  the  yield  in  kilotons,  and  P  is  the  ambient  pressure 

b 

at  the  burst  point  in  atmospheres.  H  is  the  effective  scale  height 
in  kilmeters  for  the  (upward)  direction  being  considered. 

At  distances  r  >  rQ  we  use  the  weak  shock  equations  derived  by 
Reed  (1959).  These  are: 


where  t^  is  the  positive  phase  duration  of  the  shock,  r  is  the  ratio 
of  specific  heats,  and  c  is  the  speed  of  sound.  Nonlinear  (weak  shock) 
effects  are  contained  in  Eq.  (3)  within  the  factor  t°/t+.  This  factor 
as  given  by  Eq.  (4)  includes  both  the  effects  of  dissipation  at  the  shock 
front  and  positive  phase  lengthening  due  to  the  supersonic  front  velocity. 

The  dependence  of  t+  on  values  at  may  be  eliminated  by 
taking  the  logarithmic  derivatives  of  Eqs.  (3)  and  (4)  and  comb ini  *  to 
obtain: 
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For  an  exponential  atmosphere  this  becomes: 


where 


n 


(7) 


The  value  of  n  is  obtained  from  Eq.  (1) : 


where  f'  is  the  derivative  of  f  with  respect  to  its  argument  and  n^  is 

the  value  of  n  for  a  homogeneous  atmosphere  (H  ») .  The  quantity  n^ 

is  plotted  vs.  relative  overpressure  in  Fig.  3.  Note  that,  consistent 

AP 

with  Figure  2,  the  minimum  value  of  —  corresponding  to  n  *  0  occurs  at 
r  =  that  is  at  an  altitude  3  scale  heights  above  the  burst  point. 
Combining  Eqs.  (6)  -  (8)  we  obtain: 
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At  r  =  r  ,  — —  =  .2,  and  according  to  Figure  3,  n,  =  -  — .  Thus  we 
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obtain  from  Eq.  (9): 
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FIGURE  3.  The  Parameter  n,  vs  Relative  Overpressure.  From  LEHTO  and  LARSON 


which  r -  mits  Ecs.  (3)  and  (4)  to  be  written  as: 

. .  -r  /2H  r  / r 

v  ■  ?r (a)  ro  jj +  6  °  [Ei  (Ih!  -  E±  (t§ 


r  -(r  -r)/2H 


.2  —  e  0 


1  +  e  r°/2H[Ei  (iff)  -  Ei  &)]}' 


where  rQ  is  given  by  Eq.  (2)  and  where  the  exponential  integral  is 
defined  as: 
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At  the  base  of  the  ionosphere  (considering  only  low  altitude 

-ro/2H  |i_| 

detonations)  the  term  e  ^i  I  2H  I  ^on,inate8  ttie  radicals  occurring 


in  Kas.  (11)  and  (12)  and  furthermore  may  be  approximated  as 
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Making  this  approximation  and  using  Eq.  (2)  to  eliminate 
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factors  of  rQ  e  °  we  obtain  from  Eqs.  (11)  and  (12): 
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r  /12H 

Ihe  factor  e  0  may  be  neglected  since  rQ  cannot  be  larger 

AP 

than  3H,  the  reason  being  that  =■-  ■  .2  must  occur  before  the  minimum 

“  A? 

at  r  =  3H  if  it  is  to  occur  at  all.  (According  to  Eq.  (2),  does  not 
o  *• 

become  as  sm?'l  as  .2  for  sea  level  yields  in  excess  of  60  MT.)  We  set 

Z-Zfa 

—  -  — —  >  where  z  is  the  shock  altitude*  z^  is  the  burst  altitude 

(e  "  »  P  ) ,  and  H  is  the  scale  height  in  the  vertical  direction, 

b  s 

Then  taking  Hg  =  7  km,  y  =  1.4,  and  c  *  .31  km/sec  as  representative 
values  below  120  km  altitude  we  obtain  the  final  form  of  our  expressions 
for  the  shock  positive  phase  duration  and  relative  overpressure: 


t 

’+ 


i.i  y1/3  p;1'12 

D 


*  -  .028  Y1/3  p:1/12  —A  ez/28 
P  b  r 


(16) 


(17) 


For  4MT  at  5  km  and  for  shock  propagation  at  45°  from  the  hori- 

AP 

zontal  Eq.  (17)  predicts  r~  =  1.2  at  100  km  altitude.  The  SAP  and 

*  A? 

SHELL  results  shown  in  Figure  1  both  give  about  =  1.6  at  this  alti¬ 
tude.  That  the  agreement  is  oni *  approximate  is  not  surprising  in 

view  of  the  several  assumptions  made  and  in  particular  due  to  the  fact 

AP 

that  tne  weak  shock  theory  employed  becomes  inaccurate  when  —  *  1.  For 
smaller  yields  or  propagation  closer  to  the  horizontal  the  theoretical 
results  should  be  mor^  accurate. 

The  properties  of  the  shock  when  it  reaches  the  ionospheric 
level  are  given  by  Eqs.  (16)  and  (17).  These  properties,  which  will 
be  used  in  the  subsequent  discussion,  are  as  follows: 

a.  The  relative  overpressure  and  positive  phase  duration  are 
both  proportional  to  cube  root  of  yield. 

Both  relative  overpressure  and  positive  phase  duration 
are  insensitive  to  height  of  burst  at  least  for  low  altitude 


b. 


detonations.  Taking  z  to  be  100  km  and  recalling  that 
—2  jYi 

P  o,  e  b  where  H  is  taken  to  be  7  km  we  find  that  the 
b  , 

respective  height  of  burst  dependences  are:  t+  'v  [  1  +  2^/59 j 
and  AP/P  ^  [l  +  145 J  ,  where  z^  is  in  kilometers. 

c.  For  a  given  altitude  the  positive  phase  duration  is 
independent  of  range  and  the  relative  overpressure 
decreases  linearly  with  range. 

d.  A  megaton  detonation  at  sea  level  has  a  positive  phase 
duration  of  about  39  seconds  at  100  km. 

3.2.2  Direct  Effects  of  the  Shock 

We  believe  the  direct  effects  of  the  shock  to  be  the  source 
of  the  many  observations  of  ionospheric  disturbances  following  nuclear 
tests  with  periods  between  30  seconds  and  10  minutes,  periods  of  about 
a  minute  being  predominant  (BAKER,  D.  M. ,  1968).  These  disturbances 
were  observed  by  a  doppler  technique  between  altitudes  of  150  and 
200  km  and  appeared  to  propagate  at  speeds  of  about  300  meters/sec. 
Typically  the  explosions  involved  were  the  order  of  a  megaton. 

An  example  is  shown  in  Figure  4  of  the  disturbance  observed 
following  the  United  States  megaton  range  detonation  Housatonic  on 
30  October  1962  (GLASSTONE,  S.,  1964,  pg.  677e).  BAKER  (1968)  points 
out  that  a  feature  occurs  on  4  MHz,  which  is  reflected  from  180  km, 
about  17  seconds  earlier  than  a  5  MHz,  which  is  reflected  from  190  km. 

He  finds  therefore  that  the  disturbance  v's  propagating  upward  with 
a  speed  approximately  equal  to  the  local  sound  speed.  On  this  basis 
and  because  of  the  known  acoustic  ducting  properties  of  the  atmosphere 
he  hypothesizes  that  the  acoustic  signal  has  the  bulk  of  its  energy 
confined  in  a  duct  below  110-120  km  altitude  and  that  it  is  leakage 
from  this  duct  that  is  responsible  for  the  ionospheric  disturbance. 

The  point  we  wish  to  make  is  that  the  observed  periods  of 
these  ionospheric  disturbances  are  approximately  the  same  as  the  periods 
which  the  shock  from  a  megaton  range  explosion  would  contain  on  reach¬ 
ing  the  ionospheric  level.  We  say  approximately  because  among  other 
things  it  is  not  clear  whether  the  periods  contained  in  the  shock  are 
the  order  of  2t+  or  some  slightly  larger  value.  In  terms  of  Figure  1 
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this  means  that  the  portion  of  the  shock  which  initially  propagates 
at  angles  of  from  about  20°  to  35°  is  responsible  for  these  short 
period  ionospheric  disturbances. 

Ionospheric  disturbances  with  periods  greater  than  10  minutes, 
which  we  now  discuss,  are  related  in  a  complex  manner  to  the  portions 
of  the  shock  which  initially  propagate  at  steeper  angles  to  the  hori¬ 
zontal  . 

3.2.3  Long  Period  Acoustic-Gravity  Wave  Generation  by  the  Upward 
Going  Shock 

As  previously  noted.  Figure  2  shows  a  dramatic  drop  in  relative 
overpressure  commencing  at  about  the  100  km  altitude  level.  The  reason 
for  this  can  be  seen  by  examining  the  SHELL  outputs  for  the  detailed 
flow,  one  of  which  is  shown  in  Figure  5.  This  Figure  is  a  contour  plot 
of  relative  pressure  as  a  function  of  altitude  and  horizontal  dis¬ 
tance  at  a  time  of  600  seconds.  The  spherically  expanding  and  upward 
moving  shock  wave  is  refracted  at  the  100-120  km  altitude  level  where 
the  ambient  temperature  beings  to  increase  radically.  Of  course,  the 
vertically  propagating  shock  front  is  unrefracted.  The  result  is  the 
generation  of  a  horizontally  propagating  disturbance  in  the  flow  behind 
the  shock  front.  Thi3  disturbance  drains  energy  from  the  shock  which 
continues  outward  and  upward. 

The  essential  feature  of  the  Greene-Whi taker  calculation  is 
that  the  effective  source  of  the  long  period  gravity  wave  disturbance 
is  located  at  an  altitude  of  100-120  km  above  the  burst.  The 
effective  source  commences  10  minutes  or  more  after  a  sea  level  detona¬ 
tion,  this  being  the  time  required  for  the  shock  wave  to  reach  the 
ionospheric  level  and  for  the  rotating  flow  of  Figure  3  to  become 
established.  This  delay  time  correction  to  the  travel  time  has  been 
observed  following  several  low  altitude  explosions  (HERRON,  T.J.,  1971). 

The  velocity  vector  plots,  which  are  not  shown,  indicate  that 
the  dominant  motion  in  the  disturbance  depicted  in  Figure  5  is  a  rotacing 
or  vortex  f3ow  and  on  this  basis  we  proceed  to  estimate  the  disturbance 
kinetic  energy.  The  geometry  assumed  for  this  vortex  flow  is  indicated 
in  Figure  6. 
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The  kinetic  energy  of  the  motion  is: 


\  p  v  dV, 


where  the  volume  element  for  the  region  containing  the  vortical 
flow  is: 

dV  *  2  tt  R  (2irr’  dr’) 

For  the  assumed  flow  we  may  define  a  constant  vorticity 
within  the  rotating  region: 


o)  -  h  |V  x  v  |  «  —7-“ 


in  terms  of  which  Eq.  (18)  becomes: 

Eg  =■  8  tt2u2  R  j  p(r')  r*3  dr', 


and  approximating  p(r')  by  p  ,  the  ambient  density,  this  results  in: 
2  2  4 

E  -  2  ir  u  r  R  p  (22] 

s  a 


We  estimate  R  *  150  km  and  r  -  30  km  based  on  Figure  5. 

This  latter  value  is  approximately  the  gradient  scale  of  the  tempera¬ 
ture  change  responsible  for  the  shock  refraction  and  development  of 
the  vortical  motion,  i.e.  T/VT  =  r.  We  estimate  the  maximum  flow 

velocity  v(r)  as  being  approximately  equal  to  the  sound  speed  at  130  km, 

-9 

that  is  about  .64  km/sec.  Then  for  an  ambient  density  of  7.6  x  10 
3 

kg/meter  we  obtain  E  -  5  KT,  that  is  about  .1%  of  the  total  4  MT  yield. 

For  smaller  yields  than  our  4  MT  example,  the  yield  scaling 
of  kinetic  energy  in  the  horizontally  propagating  disturbance  is  easily 
obtained.  The  kinetic  energy  is  *5  pv  V,  where  p  is  the  density,  v  the 
fluid  velocity,  and  V  the  volume  over  which  the  rotational  motion  takes 
place.  For  a  weak  shock  we  expect  only  the  fluid  velocity  to  be  a 
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function  of  the  detonation  parameters.  Since  this  velocity  is  pro- 

protional  to  the  relative  overpressure  we  find  from  Eq.  (17)  that  the 

energy  in  the  horizontally  propagating  disturbance  is  proportional  to 
2/3 

Y  .  Thus  with  a  proportionality  constant  determined  by  the  AMT  example 
we  may  write: 

E  /v  -  2  x  10"2  Yl/3.  (23) 

s/Y 

We  will  compare  this  estimate  with  the  energy  available  in  the  rising 
fireball  for  wave  generation,  which  is  derived  below. 

3.3  FIREBALL  GENERATION  OF  ACOUSTIC-GRAVITY  WAVES 

The  fireball  is  a  buoyant  rising  fluid  element  and  as  such 
deforms  into  a  torus  or  vortex  ring  configuration.  The  flow  is  illustrat¬ 
ed  in  Figure  7.  Duri* g  the  initial  rise  period,  that  is,  well  before  the 
stabilization  time  the  flow  around  this  vortex  ring  can  be  treated  as  the 
flow  around  an  equivalent  solid  object.  Warren  (1960)  has  calculated 
the  wave  emission  and  wave  drag  on  a  sphere  moving  vertically  in  a 
stratified  incompressible  fluid.  His  theory  predicts  that  gravity  wave 
emission  is  negligible  when 


where  ft  =  v/2a,  v  is  the  sphere  velocity,  and  a  the  sphere  radius.  The 
Brunt-Vdisdia  frequency  N,  is  the  high  frequency  cutoff  for  gravity  wave 
propagation.  The  onset  of  wave  emission  at  ft/N  =  1  has  been  experimenta- 
ly  verified  (MOWBRAY,  D.  E.  and  B.  S.  H.  Rarity,  1967). 

A  buoyant  sphere  rises  at  a  velocity  (SCORER,  R.  D.,  1950) 


v  =  (g  B  a)"*,  (25) 

where  g  is  the  acceleration  of  gravity  and  8  *  is  the  buoyancy 

defined  as  the  difference  in  density  between  the  fireball  and  its  sur¬ 
roundings  divided  by  the  ambient  density. 
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FIGURE  7.  Cutaway  of  Nuclear  Cloud  Showing 
Internal  Flow.  From  GLASSTONE 
(1964,  p.  31). 


Warren’s  criterion  for  significant  gravity  wave  emission  there¬ 
fore  gives: 


a  >  g  $ 


2 


(26) 


=  5.6  6  km , 


for  a  Brunt-Vaisala  period  of  5  minutes.  Because  of  the  very  high 
fireball  temperatures  involved  at  early  times  we  may  approximate  B  by 
unity.  Then  Eq.  (26)  states  that  a  fireball  radius  of  5.6  km  is  required 
for  the  onset  of  appreciable  gravity  wave  emission.  According  to  the 
data  in  Glasstone  (1964, p  10)  this  would  require  a  yield  in  excess  of  30MT. 
We  conclude  that  for  sea  level  detonations  only  very  large  yields  could 
result  in  gravity  wave  emission  during  the  initial  period  of  the  rise. 

When  the  fireball  reaches  its  stabilization  altitude  and 
proceeds  toward  hydrodynamic  equilibrium  with  its  surroundings,  conditions 
more  favorable  for  gravity  wave  emission  do  occur.  However,  before  con¬ 
sidering  the  details  of  how  this  wave  emission  arises  we  conclude  our 
discussion  of  the  initial  rise  period  by  computing  the  fraction  of  the 
total  yield  which  is  contained  in  the  kinetic  energy  of  the  rising  fire¬ 
ball  for  a  sea  level  detonation.  We  wish  to  compare  this  energy  with 
our  previous  estimate  for  the  energy  available  in  the  shock  for  long 
period  gravity  wave  generation.  The  fireball  kinetic  energy  is: 

Ef  *=  h  Mv2  =  h  M  g  B  a  (27) 

where  we  have  used  Eq.  (25)  for  V,  The  mass  occurring  in  Eq.  (27)  is  not 
the  mass  of  the  fireball  itself  but  rather  that  of  the  considerably 
larger  mass  of  ambient  air  which  moves  with  the  fireball,  the  so  called 
virtual  mass  (DARWIN,  1953).  This  mass  is  equal  to  half  the  mass  of 
ambient  air  which  would  occupy  the  fireball  volume.  Accordingly  we 
write: 

Ef  *  |  a4  p  g  B  (28) 
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According  to  Glass  tone  (1964,  pp.  74  -  75)  the  terminal  fireball  radius 

is  about  750  feet  for  a  20  KT  detonation.  The  temperature  at  the  time 

(1  sec)  this  radius  is  achieved  is  about  10  times  ambient,  and  hence 

assuming  pressure  equilibrium  6  ■  .9.  Assuming  cube  root  scaling  with 

3  1/3 

yield  we  may  then  write,  a  ■  8.45  x  10  Y  cm,  where  Y  is  in  kilotons. 

Eq.  (28 )  then  gives  the  fraction  of  the  yield  which  is  available  as 
kinetic  energy  of  fireball  rise: 

~  =  1.5  x  10~4  Y1/3  (29) 

Comparing  this  with  Eq.  (23)  we  find  that  for  a  yield  of  about  1.6  MT 
=  Es*  Thus  for  yields  in  excess  of  this  value  we  expect  the  energy 
available  in  the  rising  fireball  for  long  period  gravity  wave  generation 
to  be  greater  than  the  similar  energy  available  in  the  upward  going  shock. 

We  now  turn  to  a  description  of  how  this  fireball  rise  energy 
becomes  available  for  the  generation  of  long  period  gravity  waves. 

3.3.1  Oscillating  Vortex  Model  of  Fireball  Stabilization 
In  this  section,  the  relation  between  a  fireball  in  the 
terminal  phase  of  its  motion  and  the  fluid  oscillations  that  result  as 
it  stops  rising  are  considered.  An  approximate  expression  for  the 
frequency  of  the  resulting  motion  is  derived  in  terms  of  the  final  motions 
of  the  fireball  and  the  atmospheric  parameters.  Our  analysis  of  gravity 
wave  generation  by  a  rising  fireball  differs  from  the  work  of  Tolstoy 
and  Lau  (1971)  in  that  it  is  based  on  the  theory  of  buoyant  vortex  rings. 

We  consider  a  ring  vortex  where  the  core  radius  R  is  much 

-  o 

smaller  than  the  major  radius  R  as  shown  in  Figure  8.  The  substance,  a 
spheroidal  mass  of  air  surrounding  the  ring,  moves  with  the  ring.  The 
difference  between  the  mass  contained  in  the  ring  and  in  the  substance 
and  an  equal  volume  of  the  atmosphere  is: 

,f  rz2(r) 

AM  =  2  tt  J  I  (p '  —  pQ)  dz  rdr,  (30) 

0  z1(r) 

where  p'(r,z)  is  the  mass  density  inside  the  vortex,  pQ(z)  is  the  ambient 
density,  and  z^Cr) »  zi(r)  define  the  boundaries  of  the  substance. 
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On  a  contour  such  as  1,  which  contains  but  does  not  Intersect 
the  core,  the  circulation  K  is  given  by: 


where  the  dot  indicates  a  time  derivative.  The  quantity  C  depends 
logarithmically  on  the  distribution  of  vorticity  within  the  ring. 

We  assume  that  the  motion  is  similar  and  that  therefore  C  is  a  con¬ 
stant  . 

The  circulation  along  contours  such  as  1  and  2  changes  dur¬ 
ing  the  rise  according  to  the  Bjerkness  Eq: 


where  p  is  the  pressure,  taken  to  be  ambient  under  the  assumption  that 
equilibrium  has  been  obtained. 

On  a  contour  such  as  2  which  does  not  contain  the  core,  K  =  0 
at  all  time 8  and  hence  K  ■  0  at  all  times.  On  any  contour  such  as  1, 
which  contains  but  does  not  Intersect  the  core,  K  must  be  the  same.  Hence 
at  any  time  K  must  be  the  same  for  all  such  contours. 

Eq.  (32)  can  be  rewritten: 


rz2(r)  r 

*-***•-*%-/  f* +/ 


zx(r) 

(PQ  -  P’) 


z2(r) 


z^r) 


since  $  «  0.  Using  the  hydrostatic  equilibrium  equation 

po  z 

Eq.  (33)  finally,  becomes: 


K  -  g 


*2(r)  ,0-p' 


z1(r) 
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which  is  nonzero  and  independent  of  r  for  r  <  R  -  Rq  and  is  zero  for 

r  >  R  +  R  •  Since  the  model  is  intended  to  apply  to  the  later  stages 
o  t 

of  fireball  rise  when  the  buoyancy  has  become  small  we  replace  p 

in  the  denominator  of  Eq.  (34)  by  p^.  We  further  assume  that  the  fire¬ 
ball  is  small  enough  that  the  integration  in  Eq.  (34)  takes  place  over 
much  less  than  an  atmospheric  scale  height  and  that  the  factor  l/pQ  may 
be  taken  outside  the  integral.  Then  using  Eq.  (34)  we  may  rewrite 
Eq.  (30)  as: 


Conservation  of  mass  gives : 


(36) 


where  V  is  the  volume  over  i/hich  the  integration  in  Eq.  (30)  is  done 
and  p"  is  the  average  density  within  this  volume.  The  term  dV^'dt 
represents  the  volume  change  due  to  adiabatic  expansion  as  the  vortex 
moves  upwards: 


(37) 


Thus  Eq.  (36)  becomes: 


_d 

dt 


V  (p"  -  pQ)  = 


dAM 

dt 


(38) 


combining  this  result  with  Eq.  (35)  for  AM  we  obtain 


r 
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dp. 


K  +  K 


dt 


inp  R2-£iL  X-zJL  i.„° 

0  -rP2  V  P„  dt 

TTil  0 


-  0 


(39) 


making  use  of  the  fact  that 


dpo  apo  •  •  K 

•g —  z  and  z  *  C  —  this  becomes: 


dt 


".id  An  pR2  +  a2  N2  K  -  0 
K  +  K  dt 


(40) 


3P. 


where  N2  =  -  . ■"■  —~ 

Y  pc  3z 


is  the  ambient  Brunt  Vaisala  frequency 


2  CV 

and  a  =  — r  is  a  constant  of  order  unity  as  discussed  below.  Eq.  (40) 
;tk 


is  most  convenientlv  solved  in  terms  of  an  auxilliary  function  Q: 


K 


V% « 


(41) 


where : 


5  +  [.V  -  l^iL  J  Q_  o 


V7  R 


(42) 


and  where  C\fp"  R)  represents  — j  (^p"  R) . 


Assuming  that  the  bracketed  term  in  Eq.  (42)  is  slowly  varying 
we  obtain  the  WKB  solution: 
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1  N  -  — - 
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Using  Eqs.  (31)  and  (41)  this  gives  the  solution  for  z: 

f  r  2  2  <V^0  R)  i 

exp  i  /  [  a  N - — -  dt 

•  1  j  1  ^0  E 


VT0*2 


f  2  2 
[  a  N  " 


_ y^;  r 

<V%'R>  ,  * 


R 


Before  discussing  the  physical  significance  of  this  equation 


it  is  useful  to  consider  the  quantity  a  - 


in  more  detail.  The 


volume  V  is  approximately  equal  ttR  Az,  where  Az  is  the  vertical  distance 
over  which  the  z  integration  in  Eq.  (30)  is  done.  Alternatively  Az  may 
be  interpreted  as  the  effective  distance  beneath  the  ring  from  which 
fluid  is  drawi  up  to  the  ring.  Observations  of  laboratory  releases  of 
buoyant  fluid  with  a  spheroidal  substance  suggest  the  va^ue  Az  =  3R 
(TURNER,  J.  ?. ,  1960).  The  quantity  C  is  about  .13  for  a  spherical  vortex 
(LAMB,  1932).  This  gives  a  -  .62.  In  general  we  expect  a  to  be  constant 
so  long  as  C  is,  and  to  have  a  value  somewhat  less  than  unity. 

Equation  (44  )  applies  from  the  time  the  fireball  has  formed 
a  vortex  ring  and  has  entrained  sufficient  air  that  the  interior  density 
p1  is  approximately  equal  to  the  ambient  value  p  .  Thus  we  are  not  able 
to  describe  the  early  stages  of  the  rise.  However,  the  more  time  con¬ 
suming  later  stages  of  the  rise  will  occur  on  a  time  scale  the  order  of 


—  ,  that  is  on  the  order  of  the  Brunt-Vhisala  period.  For  nuclear  clouds 

which  do  not  spread  rapidly  when  approaching  their  stabilization  altitude, 

vertical  oscillations  will  occur  at  a  frequency  oN.  According  to  Eq.  (44), 

clouds  which  spread  rapidly  will  undergo  oscillations  at  a  lower  frequency 

2 

than  aN  and  the  amplitude  of  the  oscillations  will  subside  as  1/R  as  the  cloud 
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spreads.  If  the  cloud  radius  spreads  sufficiently  rapidly  (such  as  ex¬ 
ponentially  on  a  time  scale  shorter  than  )  there  will  be  no  oscillatory 
motion  and  the  approach  to  equilibrum  will  be  overdamped. 

In  any  case,  for  megaton  range  detonations  where  the  nuclear 
cloud  penetrates  into  the  very  stable  stratosphere,  a  rapid  growth  in 
cloud  radius  ensues  (GLASSTONE  1964,  pp.  35-37  ).  This  would  tend  to 
reduce  the  amplitude  of  gravity  waves  due  to  vertical  oscillations  of 
the  cloud.  In  general,  for  given  atmospheric  conditions  there  may  be 
an  optimum  yield-height  of  burst  combination  for  gravity  wave  generation 
by  the  stabilized  fireball  but  it  is  currently  unknown. 

To  gain  further  insight  into  the  wave  pattern  generated  by 
this  type  of  buoyant  oscillation,  laboratory  experiments  have  been  per¬ 
formed  at  Mt .  Auburn,  as  described  below. 

3.3.2  Laboratory  Observations  of  Gravity  Wave  Generation  by 
Buoyantly  Rising  Masses 

A  photograph  of  the  experimental  apparatus  is  shown  in  Figure  S. 
The  experiments  were  perormed  in  a  plexiglas  tank  5'  x  3'  x  3V  deep. 

A  linearly  stratified  density  fluid  (salt  and  water  )  was  established 
in  the  tank  by  carefully  introducing  fluid  layers  of  successively  de¬ 
creasing  density.  This  stratification  very  closely  approximates  an 
exponential  medium  over  the  dimensions  of  the  experimental  region.  The 
density  gradient  was  checked  by  carefully  weighing  a  300  cc  teflon  rod 
suspended  at  varying  depths  in  the  fluid,  and  after  some  days  it  was 
found  that  the  initial  density  discontinuities  had  smoothed  out  to  form 
a  linear  density  gradient.  Using  ethyl  alcohol  as  the  buoyant,  miscible 
fluid,  a  number  of  buoyant  releases  were  made  from  a  release  mechanism 
on  the  base  surface  of  the  tank.  Visualization  of  the  wave  motions  in 
the  fluid  was  achieved  using  small  neutrally  buoyant  polystyrene  beads 
which  acted  as  excellent  scattering  centers  for  a  well  collimated  light 
beam.  Thus  a  cross-section  through  the  fluid  motions  was  observed. 

A  16  mm  cine  camera  was  used  to  record  the  internal  wave  motions 
and  a  schematic  drawing  based  on  a  sample  frame  is  shown  in  Figure  10. 

For  this  release  the  buoyant  fluid  was  dyed  to  be  readily  observable. 

The  formation  of  a  vortex  ring  structure  and  vertical  oscillations  and 
rapid  spreading  of  the  initially  buoyant  fluid  as  it  reached  its  equi¬ 
librium  level  were  clearly  observed. 
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Typically  the  time  required  for  the  top  of  the  released  fluid 
to  reach  its  maximum  altitude  was  about  8  seconds,  somewhat  shorter  than 
the  Brunt-VSishlM  period  in  the  tank  (9.8  seconds). 

The  wave  pattern,  which  was  more  complex  than ^anticipated,  is 
conveniently  divided  into  the  four  regions  shown  in  Figure  10:  (1)  Directly 
above  and  below  the  source  region  small  amplitude  oscillations  of  the 
beads  were  observed.  The  period  of  these  oscillations  as  nearly  as  could 
be  determined  was  equal  to  the  Brunt-Vaisala  period.  (2  )  The  amplitude 
and  the  period  of  the  oscillations  increased  with  increasing  angle  from 
the  vertical.  The  amplitude  reached  a  maximum  at  an  angle  of  about  60° 
where  the  period  was  about  12  seconds.  (3)  The  amplitude  then  decreased 
slightly  with  increasing  angle.  This  region  existed  over  a  relatively 
small  angle  and  it  was  difficult  to  observe  changes  in  the  period  with¬ 
in  this  region.  (4)  Finally,  at  very  small  angles  to  the  horizontal  planes 
of  stratification,  motion  of  the  beads  was  observed  which  was  not  periodic, 
rather  the  beads  were  displaced  away  from  the  source  region  and  did  not 
return.  This  motion  is  clearly  related  to  the  horizontal  spreading  of  the 
initially  buoyant  fluid  as  it  equilibrates. 

The  propagation  theory  for  acoustic-gravity  waves  predicts  that 
waves  of  period  x  whose  vertical  wavelength  is  less  than  a  scale  height 
should  propagate  at  an  angle  0  from  the  vertical  (HINES,  C.O.,  1960), 
where : 


BV  , 

-  ■  COS  0 

T  1  1 


This  propagation  law  does  appear  to  be  obeyed  in  the  experiment  since 

2 

t  =  at  6  «  0,  i  and  x  ~-yyr  TBV  at  0  =  tt/3 ,  2u/3.  In  theory  there 

should  be  no  propagating  solution  exactly  at  0  *  0,  tt  where  t  «  x  .  However, 

BV  1 

because  the  characteristic  wavelengths  in  the  experiment  are  much  larger 
than  the  dimensions  of  the  tank  we  are  unable  to  distinguish  propagating 
from  spatially  damped  disturbances. 

Perhaps  the  most  interesting  aspect  of  the  laboratory  observa¬ 
tions  is  the  confirmation  of  the  theoretical  prediction  that  much  of  the 
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buoyant  rise  energy  is  utilized  to  generate  waves  in  the  period  range 
1-2  Our  expectation  is  that  an  analogous  process  of  wave  gen- 

eration  would  occur  for  a  stabilizing  nuclear  cloud.  Ibwever,  while 
we  are  able  to  estimate  the  period  of  the  largest  amplitude  portion 
of  the  disturbance  we  have  not  estimated  what  fraction  of  the  avail¬ 
able  energy  goes  into  very  long  periods,  that  is  periods  many  times 
the  Brunt-Vdisala  period.  These  long  period  waves  would  propagate 
at  shallow  angles  to  the  horizontal  and  it  may  be  that  atmospheric 
refraction  would  prevent  their  reaching  ionospheric  levels.  In  this 
connection  it  is  interesting  to  note  that  ionospheric  disturbances  of 
about  3  minutes  period,  which  are  associated  with  severe  weather,  are 
not  normally  observed  more  than  250  km  horizontally  from  the  apparent 
source  location  (GEORGES,  T.  M. ,  1968  ).  Whether  the  absence  of  longer 
periods  at  greater  horizontal  distances  is  a  result  of  properties  of 
the  source  or  of  atmospheric  filtering  and  refraction  is  not  known. 

3.4  CONCLUSIONS 

In  the  preceding  sections  we  have  attempted  to  relate  the 
different  kinds  of  disturbance  seen  in  the  ionosphere  to  the  hydrodynamic 
motions  accompanying  a  low  altitude  nuclear  explosion.  The  short  period 
ionospheric  disturbances  seem  to  be  caused  by  the  portion  of  the  shock 
front  which  is  reflected  from  the  100-120  km  altitude  level.  This 
identification  is  made  because  the  periods  contained  within  the  shock 
when  it  reaches  this  altitude  level  are  the  order  of  the  observed 
ionospheric  periods,  that  is,  about  a  minute.  We  expect  the  ionospheric 
periods  to  scale  with  explosion  parameters  in  the  same  way  as  does  the 
shock  positive  phase  duration  at  the  100-120  km  altitude  level.  The 
ionospheric  periods  should  therefore  be  proportional  to  the  cube  root 
of  yield  and  essentially  independent  of  height  of  burst. 

The  origin  of  ionospheric  disturbances  with  periods  exceeding 
10  minutes  is  a  more  complex  problem  since  both  the  upward  going  shock 
and  the  rising  fireball  are  possible  hydrodynamic  sources.  For  yields 
in  excess  of  about  a  megaton  the  energy  available  in  fireball  rise  for 
the  production  or  long  period  disturbances  exceeds  that  available  in  the 
upward  going  shock. 
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The  mechanism  described  by  Greene  and  Whitaker  (1968)  for 
the  generation  of  long  period  ionospheric  waves  involves  a  refraction 
of  the  upward  going  shock  at  the  100-120  km  altitude  level.  This 
results  in  a  transfer  of  energy  from  the  shock  to  a  horizontally  pro¬ 
pagating  disturbance  which  is  left  behind.  The  Greene-Whitaker  mechanism 
is  generally  operative  since  conditions  for  this  shock  refraction  are  a 
permanent  feature  of  the  atmosphere. 

The  fireball  is  most  efficient  in  the  generation  of  long 
period  acoustic-gravity  waves  when  it  has  reached  its  terminal  altitude 
and  is  approaching  hydrodynamic  equilibrium  with  the  atmosphere.  It 
may  then  undergo  vertical  oscillations  with  a  period  somewhat  larger  than 
the  Brunt-Vaisala  period.  The  efficiency  of  wave  generation  by  this 
mechanism  is  very  sensitive  to  ambient  conditions  at  the  tropopause  as 
well  as  to  burst  parameters.  In  particular  nuclear  clouds  which  under¬ 
go  rapid  spreading  upon  stabilization  appear  to  be  inefficient  wave 
generators.  In  theory,  and  this  is  supported  by  laboratory  experiments, 
the  terminal  phases  of  fireball  rise,  when  the  mechanism  is  efficient, 
produce  a  wave  spectrum  with  amplitudes  peaked  at  one  i  two  times  the 
Brunt-Vdisald  period.  These  periods  propagate  at  relat-ively  steep  angles 
to  the  horizontal  and  would  primarily  affect  the  ionosphere  near  the  burst. 
Whether  significant  wave  amplitudes  could  be  produced  at  much  longer 
periods  is  presently  unknown  as  is  the  question  of  whether  these  long 
periods,  which  propagate  at  shallow  angles  to  the  horizontal,  would  be 
prevented  by  atmospheric  refraction  from  reaching  ionospheric  levels. 
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4.  EXPERIMENTAL  INVESTIGATION  OF  INTERNAL  WAVES 
GENERATED  BY  A  BUOYANTLY  RISING  FLUID 


4.1  INTRODUCTION 

There  has  been  considerable  interest  in  the  generation  and 
propagation  of  waves  in  the  atmospheres  of  earth  and  other  planets,  and 
a  number  of  authors  (LAMB,  1945;  ECKART,  1960;  TOLSTOY,  1963;  PIERCE, 
1966)  have  discussed  the  mathematical  formulation  of  the  relevant  wave 
equations.  Experimental  investigations  have  also  been  carried  by 
MOWBRAY  and  RARITY,  1967;  WU,  1969;  HURLEY,  1969;  SCHOOLEY  and  HUGHES, 
1972. 

The  present  work  is  the  first  experimental  investigation  of 
the  internal  gravity  wavefield  generated  by  the  motion  of  a  buoyantly 
rising  fluid  in  an  incompressible,  density-stratified  atmosphere.  Also 
for  the  first  time  the  trajectories  of  specific  fluid  elements  in  the 
wavefield  have  been  observed  and  the  results  appear  to  be  in  close 
accord  with  theoretical  predictions. 

4.2  EXPERIMENTAL  FACILITY  AND  TECHNIQUE 
4.2.1  Description  of  Apparatus 

Figure  1  shows  a  photograph  of  the  main  experimental  facility, 
a  plexiglas  tank  of  dimensions  5'  x  3'  x  3V  deep.  This  tank  contains 
the  stratified  fluid,  which  in  our  experiments  was  a  solution  of  vary¬ 
ing  amounts  of  salt  in  water.  Distilled  water  was  used  to  improve  the 
light  transmission  characteristics  of  the  solution.  The  procedure  for 
filling  the  tank  with  the  density-stratified  solution  was  as  follows: 
first  a  solution  with  the  maximum  density  of  salt  was  filled  into  the 
tank,  typically  to  a  depth  of  3  or  4  inches.  Next  a  similar  volume  of 
somewhat  lower  density  solution  was  floated  on  top  of  the  first.  This 
was  achieved  by  flowing  the  new  solution  onto  a. thin  board  floating  on 
the  top  surface  of  the  iquid  already  in  the  tank  so  as  to  minimize 
vertical  mixing.  Using  this  technique,  adjacent  fluid  strata  with 
density  differences  of  less  than  5x10  could  be  clearly  established. 
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The  operation  was  then  repeated  with  a  succession  of  fluid  strata  until 

the  desired  overall  density  deficit  from  top  to  bottom  was  obtained. 

After  a  period  of  some  days  ionic  diffusion  at  the  interfaces 

between  successive  strata  smoothed  out  the  density  discontinuities  and 

eventually  provided  a  linear  density  profile  in  the  vertical  direction. 

The  density  profile  was  measured  by  carefully  weighing  a  teflon  cylinder 

(300  cc)  submersed  at  varying  depths  in  the  stratified  fluid.  A  density 

change  Ap  in  the  local  average  value  p  can  be  determined  to  an  accuracy 
Ap  AW  ° 

-±L  *=  rp  where  W  is  the  apparent  weight  of  the  teflon  cylinder  suspended 
p  w 
o  o 

in  the  fluid.  This  allowed  the  density  to  be  determined  t-  b-.t-L  r  than 
lxl0~4. 

Experiments  were  conducted  in  a  number  of  test  .ids  whose 


density  stratifications 
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—  r —  / ,  ranged  between  5x10  cm  and 
PQ  J 
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5x10  cm  .  Thus  the  Brunt-Vaisala  periods,  =  2ir/(ag)  ,  for  the  strati¬ 
fied  fluid  mixes  used  in  our  experiments  were  in  the  range  between  3 
and  10  seconds. 

Internal-gravity  waves  were  excited  in  the  stratified  test 
mix  by  one  of  two  methods.  In  a  series  of  preliminary  experiments  a 
sphere  was  suspended  in  the  test  tank  and  using  a  small  D.C.  motor  the 
sphere  was  allowed  to  oscillate  vertically  at  a  controlled  amplitude 
and  frequency.  Measurements  were  made  of  the  amplitude  and  frequency 
of  waves  as  a  function  of  angle  from  the  axis  of  motion  of  the  oscillat¬ 
ing  sphere.  The  results  of  these  experiments  are  reported  in  Section  4,4, 
below. 

The  main  series  of  experiments  in  this  investigation  were 
conducted  by  releasing  a  small  volume  of  buoyant  fluid  from  a  control 
mechanism  on  the  floor  of  the  plexiglas  tank.  The  release  mechanism 
has  a  facility  for  refilling  with  buoyant  fluid  from  beneath  the  tank; 
the  buoyant  fluid  is  then  retained  in  a  hemispherical  volume  by  a  slid¬ 
ing  lid  on  top.  The  lid  is  opened  using  a  draw  cord  which  is  led 
through  guides  on  the  tank  floor  and  up  the  end  walls  of  the  tank  to 
allow  manual  control  from  outside  the  tank.  The  lid  is  spring  loaded 
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to  allow  a  controlled  opening  time.  Ethyl  or  methyl  alcohol  were 
generally  used  as  the  buoyant  fluid  due  to  their  low  density  and  good 
miscibility  with  water,  and  in  some  cases  the  alcohol  was  dyed  to 
enhance  visualization  of  the  source  region. 

4.2.2  Internal-Gravity  Wave  Visualization 

The  wavefields  generated  in  the  test  tank  were  visualized  by 
using  neutrally  buoyant  polystyrene  beads  dispersed  through  the  strati¬ 
fied  fluid  mixture.  These  beads  were  purchased  in  the  unexpanded  state 
but  were  then  thermally  treated  to  produce  varying  degrees  of  expansion 
in  accordance  with  the  requirements  for  a  particular  experiment.  Samples 
of  beads  within  quite  narrow  density  ranges  are  thus  obtained  and  an 
even  distribution  of  beads  can  be  dispersed  throughout  the  density 
stratified  mixture.  Wave  propagation  in  the  stratified  test  fluid  may 
then  be  observed  by  examining  the  motion  of  individual  beads  which 
characterize  the  motion  of  the  surrounding  fluid. 

The  wavefield  created  is  a  three-dimensional  one  and  in  order 
to  simplify  the  interpretation  cf  data,  observations  were  made  on  a 
thin  two-dimensional  characteristic  section  of  the  wavefield.  This 
was  done  by  conducting  the  experiments  in  a  darkened  laboratory  with 
a  collimated  light  beam  illuminating  only  a  thin  vertical  section  of 
the  test  fluid  through  the  axis  of  symmetry  for  the  experiment.  The 
polystyrene  beads  act  as  excellent  light  scattering  centers  and  their 
movements  were  recorded  using  a  16  mm  cine  camera. 

4 . 3  THEORY 

4.3.1  Fundamental  Equations 

Internal  gravity  waves  in  an  incompressible  stratified  fluid 
are  governed  by  hydrodynamic  equations: 

Dun- 
P  ^  =  -  Vp  -  gpez 
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Let  p  *  pq(z)  +  P^»  P  *  Pg(z)  +  where  p^,  p^  are  small. 
Then  for  small  displacements 
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where  S  is  the  particle  displacement  from  its  undisturbed  position. 
From  (3)  and  (4)  it  follows  that  to  first  order. 


PX  +  5  •  Voo  =  0 


From  (1)  it  similarly  follows  that 
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p^  may  now  be  eliminated  from  (5)  and  (6)  so  that 
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From  Eqs.  (11)  it  follows  that' 
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and  that  o)  and  k  are  related  by  the  dispersion  relation 
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For  simplicity  we  will  consider  and  to  be  zero.  Then 
the  dispersion  relation  (13)  may  be  written  as 
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Figure  2 
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Sketch  of  group  velocity  C  versus  frequency 

o 

for  a  fixed  group  velocity  direction 
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where  A  and  <f>o  are  constant  for  any  given  fluid  particle. 
Writing 


(19a) 


(19b) 


(19c) 


we  derive  from  Eqs.  (14),  (15),  and  (18)  that  the  particle  path  is 
given  by 


,x2  ,  2.  r2  ,  2.2  .  ,  _  2.2 

(6  +  Y  )  Cx  +  a  ?2  ~  2ayZxZz  *  a  6 


(20) 


which  is  the  equation  of  an  ellipse.  As  illustrated  in  Fig.  3a,  the 
direction  of  the  major  axis  of  this  ellipse  is  in  general  not  coinci¬ 
dent  with  the  group  velocity  vector,  C  .  The  inclination  $  of  the 
major  axis  of  the  ellipse  with  the  vertical  may  be  defined  by 
(0  <  $  <  tt/2) 
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(a)  General  case;  (b) 


(21) 
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and  the  ratio  L^/L^,  w^ere  L2»  are  semi-major  and  semi-minor 
axes  of  the  ellipse,  respectively,  may  be  written  as 


2  2  2  2  2  2 

(cT  +  y  +  6  )  cos  2j>  +  (6  +  y  -  a) 
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Substituting  from  (19)  into  (21)  we  find  that 
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Theoretical  arguments  not  presented  here  indicate  that  a 

source  of  frequency  w  radiates  most  strongly  near  the  angle 

e  =  ±  cos-i(w/u)  )  (PIERCE,  A.  D.,  1963).  At  precisely  this  angle, 
o 

Eq.  (17)  shows  the  group  velocity  is  zero.  However,  at  angles  only 

slightly  different  from  this  value,  the  group  velocity  has  a  finite 

2  2  2 

value  and  the  amplitude  is  expected  to  be  maximum.  For  w  -  u)q  cos  6, 
Eq.  (23)  reduces  to: 


tan  2<j>  *  tan  26 


which  is  satisfied  by 
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Equation  (22)  may  be  Simplified  and  in  general  can  be  written 
as 


2  2  2 

Thus  in  the  case  where  co  -*■  co  cos  6,  then  L,  -*•  0  and  the 

o  1 

particle  trajectory  becomes  a  straight  line  whose  orientation  is  in 

the  direction  of  0  (=  <p)  .  See  Fig.  3b. 

In  the  special  case  of  0  =  0,  the  only  propagating  frequency 

is  coq,  so  Lj  0,  and  $  =  0°.  Since  the  straightline  orbits  are  in 

the  direction  of  <[>  they  will  thus  be  aligned  in  the  vertical  direction. 

See  Fig.  3c.  At  the  other  extreme  case  of  horizontal  propagation,  as 

0  -♦>  90°,  <p  may  go  either  0°  or  90°  depending  on  whether  co^  is  less  than 
2 

or  greater  than  <oq/2  and  elliptical  trajectories  should  in  general  be 
observed  in  the  horizontal  direction. 

2  2 

Combining  Eqs.  (23)  and  (26)  we  can  eliminate  co  /to  to  obtain 


Thus  if  one  measures  $  and  the  ratio  L^/Lj  for  an  observed  particle 

trajectory,  it  is  possible  to  compute  0  without  explicitly  taking  the 

frequency  into  account.  (This  presumes  0  <  0  <  90°,  0  <  <p  <  90°.). 

One  other  result  which  should  be  mentioned  is  the  sense  in 

which  particles  move  around  elliptical  orbits.  One  may  note  from 

Eqs.  (18)  that,  given  co^  <  co^,  k  >  0,  k  <  0  (C  >  0,  C  >0)  that 

o  x  z  gx  gz 

d£z/dt  is  negative  when  has  its  maximum  value.  In  this  case  the 
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orbits  would  appear  clockwise  if  the  source  is  to  the  lower  left  of 
the  observer.  In  a  similar  manner,  the  orbits  would  be  counter  clock¬ 
wise  if  the  source  is  to  the  lower  right. 

4.4  RESULTS 

4.4.1  Bobbing  Sphere  Experiments 

A  series  of  preliminary  experiments  was  carried  out  employ¬ 
ing  an  oscillating  sphere  as  the  localized  disturbance  in  a  stratified 
fluid.  A  schematic  illustration  is  presented  in  Fig.  4.  Sphere  sizes 
from  0.5-  to  2.5-inch  diameters  were  used  and  oscillation  frequencies 
from  approximately  one-third  to  one  times  the  Brunt-Vaisala  frequency 
were  examined.  A  16  mm  cine  camera  recorded  the  movement  of  the  sphere 
and  the  wave  patterns  indicated  by  the  particle  motions,  viewing  only  a 
vertical  two-dimensional  plane  as  described  in  Section  4.2. 

Analysis  of  the  data  consisted  of  plotting  the  movement  of 
the  neutrally  buoyant  beads  at  various  locations  in  the  wavefield. 

The  wave  pattern  was  clearly  observed  to  take  the  form  of  symmetrical 
"jet"  regions  on  either  side  of  the  oscillating  source,  in  the  shape 
of  an  X  as  reported  by  MOWBRAY  and  RARITY,  1967.  The  angle  of  inclina¬ 
tion  of  the  gravity  wave  jet  regions  was  found  to  vary  closely  as 

_  t 

6  =  cos  "(ea/w  )  for  values  of  u  <  id  where  0  is  the  jet  inclination 
o  o 

angle  to  the  vertical.  This  is  consistent  with  Hurley's  results 
(HURLEY,  D.  G.,  1969)  based  on  the  theory  of  an  oscillating  cylinder 
in  a  stratified  fluid.  The  small  change  in  the  Brunt-Vaisala  frequency 
through  the  density  stratified  medium  is  neglected  here  and  a  value  of 
corresponding  to  its  median  value  is  used.  Table  1  presents  data 
from  some  of  the  measurements  showing  the  correlation  between  0  and 
cos  (u)/wo)  for  a  stratified  fluid  whose  Brunt-Vaisala  period  was 
5.75  secs. 

A  detailed  analysis  was  carried  out  on  the  data  from  one  of 
our  oscillating  sphere  experiments  where  the  actual  trajectories  of  a 
number  of  the  neutrally  buoyant  beads  were  plotted  using  a  stop  frame- 
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Figure  4.  Schematic  for  Bobbing  Sphere  Experiment 


TABLE  1 


CORRELATION  OF  MEASURED  GRAVITY  WAVE  JET  ANGLE  WITH 
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analyst  projector.  All  the  particle  motions  were  found  to  follow 
straightline  paths,  and  Fig.  5  shows  displacement/time  plots  for  two 
such  particles  in  the  gravity  wave  jet,  spaced  azimuthally  but  at 
equal  radial  distances  from  the  oscillating  source.  These  clearly 
indicate  the  phase  change  across  the  jet  dimension  and  the  aggrega¬ 
tion  of  such  data  points  across  the  jet  at  a  fixed  radial  distance  from 
the  source  gives  a  measurement  of  the  wavelength.  In  our  experiments 
this  wavelength  was  found  to  agree  within  15%  with  the  diameter  of  the 
oscillating  source  region  which  was  defined  as  the  sphere  diameter 
plus  the  amplitude  of  the  imposed  oscillation. 

The  displacement/time  plots  for  the  individual  particles 
give  measurements  of  the  wave  period  and  amplitude  at  all  observed 
parts  of  the  jet  region.  The  wave  amplitude  tended  to  be  maximum  at 
the  center  region  of  the  jet,  falling  off  fairly  slowly  at  first  away 
from  the  jet  center  line  and  then  quite  rapidly  at  the  jet  edges.  The 
decay  of  wave  amplitude  with  increasing  range  R  from  the  source  seems 
to  fall  off  approximately  as  1/Jr  as  might  be  anticipated  on  the 
grounds  of  the  following  simple  consideration:  it  is  observed  that  the 
gravity  wave  jet  thickness  L  and  angle  of  inclination  0  to  the  vertical 
are  constant.  Thus  from  Eq.  (17)  the  group  velocitv  is  constant. 

Energy  conservation  requires  EC  S  =  const.,  where  F  the  energy  density 

s 

in  the  jet  and  S  is  the  area  of  the  jet  normal  to  -..e  group  velocity. 

S  =  2uRL  sin  0  and  the  energy  density  will  vary  as  A4-  where  A  is  the 
gravity  wave  amplitude.  Thus  one  would  expect  A  to  vary  as  1/^R  . 

Data  are  presented  in  Fig.  6  showing  measurements  of  the  wave  ampli¬ 
tude  in  the  jet  at  different  ranges  R  from  the  bobbing  sphere.  Also 

_j r 

shown  is  the  predicted  variation  as  R 

The  measured  periods  for  all  particles  examined  fell 
within  the  range  11.2  -  12.1  seconds  for  a  mean  vaiue  of  11.6  seconds. 

The  known  source  period  in  this  case  was  11.5  seconds.  The  average 
angle  of  inclination  of  the  gravity  wave  jet  region,  defined  by  the 
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Figure  5.  Displacement/time  plots  for  two  particles  in  the  gravity  wave  jet  in 

bobbing  sphere  experiment.  Both  particles  are  at  same  radial  distance 
from  sphere  but  are  separated  azimuthally.  Note  the  phase  difference. 


bounds  within  which  measurable  particle  displacements  were  observed, 

was  58°  from  the  vertical.  This  agrees  with  the  value  of  58°  derived 

from  cos  m/io  for  the  same  experiment, 
o 

It  was  also  found  that  the  inclination  of  the  straightline 
particle  paths  likewise  averaged  58°.  This  is  consistent  with  the 
theoretical  analysis  of  Section  4.3  which  predicts  straightline  tra¬ 
jectories  aligned  in  the  group  velocity  direction  if  io  *  cos  8. 
Illustrative  data  are  presented  in  Table  2. 

4.4.2  Buoyant  Release  Experiments 

A  number  of  experiments  were  carried  out  to  visualize  the 
wavefield  generated  by  the  motion  of  a  buoyantly  rising  volume  of 
fluid  in  a  stratified  environment.  Most  of  the  data  analyzed  was  for 
experiments  where  a  small  test  volume  of  alcohol  was  released  at  the 
tank  bottom  into  a  stratified  salt  solution  whose  stratification,  a, 
was  5x10  ^  cm  The  same  visualization  technique  was  employed  as 
described  above  where  a  cine  film  was  taken  of  the  motion  of  neutrally 
buoyant  beads  in  a  two-dimensional  vertical  plane  passing  through  the 
release  point  of  the  buoyant  fluid. 

Due  to  the  rapid  entrainment  by  the  buoyant  fluid  of  the 
heavier  brine  solution,  stabilization  of  the  rising  fluid  tended  to 
occur  at  relatively  low  levels  In  the  test  tank.  Thus  observations 
were  mainly  made  on  the  upward  directed  waves  which  provided  a  longer 
wave  path  before  reflections  at  the  tank  side  walls  or  the  top  surface 
of  the  stratified  fluid  interfered  with  the  observations. 

Whereas  in  the  experiments  with  the  bobbing  sphere  when  an 
observer  could  clearly  distinguish  the  wave  "jet"  regions  of  the  wave- 
field,  in  the  case  of  a  buoyant  fluid  release  the  wave  pattern  appeared 
considerably  mere  complex.  General  observations  were  made  on  a  number 
of  buoyant  releases  and  some  were  analyzed  in  detail;  results  are  pre¬ 
sented  below. 
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TABLE  2 


MEASURED  VALUES  OF  THE  INCLINATION  OF  PARTICLE  STRAIGHTLINE 
TRAJECTORIES  IN  A  BOBBING  SPHERE  EXPERIMENT.  cos~1(u/uq)  =  58° 


Particle 

Number 

Measured 
Inclination,  0 

Particle 

Number 

Measured 
Inclination,  6 

1 

58 

9 

63 

2 

56 

10 

64 

3 

59 

11 

57 

4 

55 

12 

58 

5 

54 

13 

57 

6 

56 

14 

59 

7 

61 

15 

55 

8 

64 

_ “ _ 
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Distribution  of  Gravity  Wave  Frequency 

In  the  case  c-t  the  bobbing  sphere  experiments,  the  period 
for  all  particles  observed  within  the  wave  jet  region  appeared  tc  be 
constant  and  equal  to  the  sphere  oscillation  frequency,  within  the 
scatter  of  the  measurements.  With  the  buoyant  fluid  release,  however, 
a  wide  spectrum  of  wave  frequencies  was  observed.  No  discernible 
gravity  wave  motions  were  detected  during  the  rise  of  the  buoyant 
fluid,  the  waves  observed  seemed  to  be  generated  subsequent  to  the 
rise,  emanating  from  the  oscillatory  motion  of  the  fluid  around  its 
stabilization  height.  Particle  trajectories  were  plotted  for  various 
parts  of  the  wavefield  and  these  showed  periods  varying  from  the  Brunt- 
Vaisala  period  to  three  to  four  times  the  Brunt  period  t  v  waves  propa¬ 
gating  at  angles  close  to  the  horizontal.  It  became  difficult  to 
resolve  waves  with  periods  longer  than  this  which  may  propagate  at 
even  shallower  angles  but  there  was  evidence  of  a  non-periodic  motion 
in  the  horizontal  direction,  assumed  to  be  associated  with  the  hori¬ 
zontal  spread  of  the  buoyant  plume  at  stabilization  height. 

The  main  reason  for  the  differences  between  the  wavefields 
generated  by  the  bobbing  sphere  and  by  the  buoyant  release  is  that  the 
source  region  in  the  latter  case  represents  an  inhomogeneous  and 
transient  source.  An  a  result,  during  the  oscillatory  phase  of  the 
stabilizing  fluid  the  trajectory  of  an  individual  particle  in  the  wave- 
field  may  be  observed  to  vary  in  time.  For  example,  a  particle  execu¬ 
ting  an  elliptical  orbit  may  subsequently  change  to  a  different  ellip¬ 
tical  orbit,  or  may  for  a  time  be  observed  to  follow  a  straightline 
path.  Figure  7  shows  sketches  of  some  actual  particle  trajectories 
recorded  for  a  buoyant  release  experiment. 

The  overall  observation  time  for  most  of  our  buoyant  release 
experiments  was  usually  four  or  five  times  the  Brunt-Vaisala  period 
for  the  stratified  fluid.  Since  the  buoyantly  rising  test  fluid 
appeared  to  stabilize  after  approximately  two  oscillations  at  close  to 
Brunt  frequency,  this  meant  that  our  observation  time  lasted  typically 
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two  to  three  periods  after  the  main  source  motion  had  ceased.  Towards 
the  end  of  this  observation  time  there  were  generally  indications,  and 
in  one  experiment  very  clear  evidence,  that  the  particle  motions  even¬ 
tually  tended  to  straightline  paths  in  the  vertical  direction  at  the 
Brunt-Vaisala  frequency.  Such  a  case  is  illustrated  by  the  particle 
trajectory  in  Fig.  7d.  This  is  in  accord  with  the  theoretical  argument 
that  the  strongest  later  arrivals  from  a  broadbased  transient  source 
should  be  at  Brunt's  frequency  since  the  group  velocity  goes  to  zero 
with  zero  slope  at  that  frequency.  Also  the  vertical  trajectories  at 
to  =  wo  are  consistent  with  the  analysis  of  Section  4.3. 

The  direction  of  the  particle  motions  around  their  elliptical 
orbits  was  observed  to  be  clockwise  or  counter  clockwise  in  agreement 
with  the  predictions  of  Eqs.  (18),  as  described  in  Section  4.3. 

It  follows  from  Eq.  (26)  that  if  a  particle  trajectory  is  a 

2  2 

straightline  path  (i.e.  L,  =0),  then  w  =  u  cos  <j>.  This  is  borne 

i  c 

out  by  the  data  presented  in  Fig.  8  where  the  data  points  relate  only 
to  particles  observed  to  be  executing  straightline  paths  in  the  wave- 
field  generated  by  a  buoyant  fluid  release. 

In  the  more  general  case  of  an  elliptical  particle  trajec¬ 
tory,  such  as  shown  in  Fig.  7a,  <{>  and  L^/ may  be  measured.  Then 
Eq.  (27)  provides  the  value  of  0  whic..  indicates  the  direction  of  the 
source  region.  From  a  number  of  such  observations  the  location  of  the 
source  may  of  course  be  determined  by  the  intersection  of  the  "0" 
lines.  In  our  experiments  the  volume  dimensions  of  the  oscillating 
source  region  are  not  specifically  defined  and  it  was  not  possible  to 
critically  assess  the  accuracy  of  0  values  calculated  from  Eq.  (27), 
However,  in  general  the  predicted  values  were  in  very  good  agreement 
with  expectations  based  on  a  knowledge  of  the  general  location  of  the 
source  region.  Quantitative  correlation  was  difficult  because  our 
observation  field  only  extended  to  about  eight  or  ten  times  the  source 
region  dimensions.  Thus  a  line  joining  a  particle  location  in  the  wave- 
field  to  the  source  volume  could  encompass  a  fair  range  of  values  for  0. 
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Description  of  Gravity  Wave  Amplitude 

Due  to  the  transient  nature  of  the  gravity  wavefield  genera¬ 
ted  in  a  buoyant  release  experiment,  it  was  not  feasible  to  confirm  the 
decay  of  amplitude  with  increasing  distance  from  the  source  region,  as 
was  done  for  the  bobbing  sphere  experiments.  However,  there  is  no 
apparent  reason  why  the  decay  should  differ  markedly  from  the  1//r  form 
observed  in  the  bobbing  sphere  experiments. 

Data  was  obtained  to  indicate  the  variation  of  gravity  wave 
amplitude  with  azimuthal  angle.  Figure  9  presents  data  showing  the 
variation  of  the  vertical  component  of  gravity  wave  amplitude  with 
azimuthal  angle.  This  vertical  displacement  is  proportional  to  the 
potential  energy  in  the  wave.  It  is  seen  that  the  vertical  displace¬ 
ment  is  a  maximum  directly  above  the  oscillating  source  region  (where ^ 

however,  the  group  velocity  C  falls  to  zero)  and  falls  to  a  minimum 

S 

in  the  horizontal  direction. 

Rise  Time  to  Stabilization 

Table  3  shows  the  results  of  measurements  made  of  the  rise 
time  for  buoyant  releases.  The  rise  time  is  defined  as  the  time 
interval  from  release  of  the  buoyant  fluid  until  the  top  surface  of 
the  stabilizing  fluid  is  first  observed  to  fall  back  in  an  oscillation 
about  the  stabilization  height.  Typically  the  stabilizing  fluid  is 
observed  to  undergo  about  two  oscillations  before  its  vertical  motion 
finally  comes  to  rest. 

It  is  seen  that  the  rise  time  appears  to  be  about  0.85  times 
the  Brunt-Vaisala  period  for  the  stratified  fluid. 


TABLE  3 


RISE  TIME  TO  STABILIZATION  OF  BUOYANT  FLUID  RELEASES 


Release 

Number 

Measured  Rise 

Rise  T  (secs) 

K 

T/T 

R  o 

11-3 

8.0 

0.82 

11-4 

8.3 

0.85 

11-5 

8.3 

0.85 

12-2 

8.0 

0.82 

11-3 

8.3 

0.85 

12-5 

8.3 

0.85 

12-6 

8.0 

0.82 

12-7 

1 

9.0 

0.92 

4.5  CONCLUSIONS 

Experiments  have  been  carried  out  to  investigate  the  internal 
gravity  waves  generated  in  an  incompressible  density  stratified  medium 
(a)  by  the  forced  oscillation  o.  «oherical  body  suspended  in  the 
stratified  medium,  and  (b)  by  the  motion  of  a  buoyantly  rising  miscible 
fluid.  The  gravity  wave  motions  were  analyzed  by  observing  the  detailed 
motion  of  small  neutrally  buoyant  beads  dispersed  throughout  the 
stratified  medium,  and  agreement  was  found  between  the  observations 
and  the  theoretical  predictions  presented  in  Section  4.3.  Observa¬ 
tions  were  also  made  on  the  rise  to  stabilization  of  the  buoyant 
miscible  fluid,  and  it  was  found  that  the  rise  time  to  stabilization 
height  was  about  0.85  times  the  Brunt-Vaisala  period  for  the  strati¬ 
fied  fluid. 

The  gravity  wavefield  generated  by  the  motion  of  a  buoyant 
fluid  was  observed  to  be  inhomogeneous  and  transient  in  nature.  How¬ 
ever,  the  theory  permits  identification  of  the  location  of  an  unseen 
source  of  internal  gravity  waves,  if  measurements  are  made  of  only 
two  characteristics  of  fluid  particles'  trajectories,  viz:  the 
inclination  of  the  major  axis  of  particle  motion,  and  the  ratio  of 
major  to  minor  axis  for  its  (in  general)  elliptical  orbit. 


5.  REPORT  CONCLUSIONS 


Conclusions  are  given  below  for:  (1)  Modeling  of  nuclear 
sources  of  acoustic-gravity  waves  producing  ionospheric  disturbances. 

(2)  The  experimental  investigation  of  internal  waves  produced  by  the 
buoyant  rise  of  a  fluid  element. 

5.1  MODELING  OF  NUCLEAR  SOURCES 

(a)  Observations  of  short  period  ionospheric  disturbances 
following  nuclear  explosions  seem  to  be  caused  by  the  portion  of  the 
shock  front  which  is  reflected  from  the  100-120  km  altitude  level.  This 
identification  is  made  because  the  periods  contained  within  the  shock 
when  it  reaches  this  altitude  level  are  the  crder  of  the  observed  iono¬ 
spheric  periods,  that  is,  about  a  minute. 

(b)  Because  of  the  above  identification  we  expect  the  short 
periods  to  scale  with  explosion  parameters  in  the  same  way  as  does  the 
shock  positive  phase  duration  at  the  100-120  km  altitude  level.  The 
ionospheric  periods  should  therefore  be  proportional  to  the  cube  root 
of  yield  and  essentially  independent  of  height  of  burst. 

(c)  The  origin  of  ionospheric  disturbances  with  periods  exceed¬ 
ing  10  minutes  is  a  complex  problem  since  both  the  upward  going  shock 
and  the  rising  fireball  are  possible  hydrodynamic  sources.  The  mechanism 
described  by  GREENE  and  WHITAKER,  1968,  for  the  generation  of  long 
period  ionospheric  waves  involves  a  refraction  of  the  upward  going  shock 
at  the  100-120  km  altitude  level.  This  results  in  a  transfer  of  energy 
from  the  shock  to  a  horizontally  propagating  disturbance  which  is  left 
behind.  The  Greene -Whitaker  mechanism  is  generally  operative  since 
conditions  for  this  shock  refraction  are  a  permanent  feature  of  the 
atmosphere.  However,  for  yields  in  excess  of  about  a  megaton  the  energy 
available  in  fireball  rise  for  the  production  of  long  period  disturbances 
exceeds  that  available  in  the  upward  going  shock. 

(d)  The  fireball  is  most  efficient  in  the  generation  of  long 
period  acoustic-gravity  waves  when  it  has  reached  its  terminal  altitude 


and  is  approaching  hydrodynamic  equilibrium  with  the  atmosphere.  It  may 
then  undergo  vertical  oscillations  with  a  period  somewhat  larger  than 
the  Brunt-Vhisala  period.  The  efficiency  of  wave  generation  by  this 
mechanism  is  very  sensitive  to  ambient  conditions  at  the  tropopause  as 
well  as  to  burst  parameters.  In  particular  nuclear  clouds  which  undergo 
rapid  spreading  upon  stabilization  appear  to  be  inefficient  wave  genera¬ 
tors.  In  theory,  and  this  is  supported  by  laboratory  experiments,  the 
terminal  phases  of  fireball  rise,  when  the  mechanism  is  efficient,  pro¬ 
duce  a  wave  spectrum  with  amplitudes  peaked  at  one  to  two  times  the 
Brunt-Vaisala  period.  These  periods  propagate  at  relatively  steep  angles 
to  the  horizontal  and  would  primarily  affect  the  ionosphere  near  the  burst, 
burst.  Whether  significant  wave  amplitudes  could  be  produced  at  much 
longer  periods  is  presently  unknown  as  is  the  question  of  whether  these 
long  periods,  which  propagate  at  shallow  angles  to  the  horizontal,  would 
be  prevented  by  atmospheric  refraction  from  reaching  ionospheric  levels. 

5.2  EXPERIMENTAL  INVESTIGATION  OF  INTERNAL  WAVES 

Experiments  were  carried  out  to  simulate  the  gravity  wavefield 
generated  by  a  fireball's  rise  to  stabilization.  In  essence,  the  experi¬ 
ment  consisted  of  the  release  of  a  small  volume  of  a  buoyant  fluid  at  the 
bottom  of  a  density  stratified  tank.  The  fluid  was  observed  to  rise  to 
stabilization  height  where  it  typically  underwent  about  two  oscillations 
at  near  the  Brunt  frequency  and  gravity  waves  generated  by  the  motion 
of  the  buoyant  fluid  were  observed  using  neutrally  buoyant  beads  dispersed 
throughout  the  density  stratified  fluid. 

The  gravity  wavefield  was  found  to  be  generated  mainly  by  the 
oscillatory  motion  of  the  fluid  at  stabilization  height  and  was  observed 
to  be  inhomogeneous  and  transient  in  nature.  The  distribution  of  gravity 
wave  frequency  and  amplitude  with  azimuthal  angle  were  identified,  and 
for  comparison  similar  data  was  discussed  for  experiments  where  the  source 
was  a  set  frequency  oscillating  sphere.  The  detailed  motion  of  individual 
fluid  elements  was  observed  to  agree  with  theoretical  predictions  presented 
in  Section  4.3.  It  was  also  shown  that  if  measurements  are  made  on  only 
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two  characteristics  of  fluid  particle  trajectories,  viz:  the  inclination 
of  the  major  axis  of  the  particle  motion,  and  the  ratio  of  major  to  minor 
axis  for  its  (in  general)  elliptical  orbit,  then  the  location  of  an 
unseen  source  of  the  gravity  waves  may  be  identified. 
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